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Abstract 



The spectrum of Hamiltonian (Markov matrix) of a multi-species asymmetric simple 
exclusion process on a ring is studied. The dynamical exponent concerning the relaxation 
time is found to coincide with the one-species case. It implies that the system belongs to 
the Kardar-Parisi-Zhang or Edwards- Wilkinson universality classes depending on whether 
the hopping rate is asymmetric or symmetric, respectively. Our derivation exploits a poset 
' structure of the particle sectors, leading to a new spectral duality and inclusion relations. 

I The Bethe ansatz integrability is also demonstrated. 

1 Introduction 

In recent years, intensive studies on non-equilibrium phenomena have been undertaken through 
00 I a variety of stochastic process models of many particle systems [SI]. Typical examples are 

' driven lattice gas systems with a simple but nonlinear interaction among constituent parti- 

cles. The asymmetric simple exclusion process (ASEP) is one of the simplest driven lattice 
gas models proposed originally to describe the dynamics of ribosome along RNA [MGP]. In 
I the ASEP, each site is occupied by at most one particle. Each particle is allowed to hop to its 

■ nearest neighbor right (left) site with the rate p (q) if it is empty. The ASEP admits exact 

analyses of non-equilibrium properties by the matrix product ansatz and the Bethe ansatz. 
^ \ The matrix product form of the stationary state was first found in the open boundary case 

^ ■ [DEHP]. Similar results have been obtained in various driven lattice gas systems in one di- 

^ ■ mension with both open and periodic boundary conditions [BE]. Applications of the Bethe 

ansatz [B] to the ASEP have also been successful. See for example the works [GS, K], [S2] 
and [DE] for the studies of the ASEP under periodic, infinite and open boundary conditions, 
respectively. In the periodic boundary case, the dynamical exponent z of the relaxation time 
to the stationary state isz = 3/2iip^q and z = 2 'd p = q according to the analysis of 
the Bethe equation. This implies that the ASEP belongs to the Kardar-Parisi-Zhang (KPZ) 
or Edwards- Wilkinson (EW) universality classes depending on whether the hopping rate is 
asymmetric or symmetric, respectively. 
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In this paper we consider the following multi-species generalization of the ASEP on a ring 
of length L. A local state on each site of the ring assumes N states 1,2,...,N. A nearest 
neighbor pair of the states is interchanged with the transition rate 



We regard the local state a = 1 as a vacant site and a = 2, . . . ,N as a site occupied by 
a particle of the ath kind. The dynamics is formulated in terms of the master equation 
= H\P{t)) on the probability vector \P{t)). We call the model the (A'" — l)-species 
ASEP or simply the multi-species ASEP. The usual ASEP corresponds to A^ = 2. 

The X Markov matrix H will be called the "Hamiltonian" although it is not 
Hermitian for p ^ q. (At p = q, it is Hermitian and coincides with the sl{N) invariant 
Heisenberg Hamiltonian.) Its eigenvalue with the largest real part is corresponding to the 
stationary state. The other eigenvalues contribute to the relaxation behavior through \P{t)) = 
e*^|P(0)). Especially those with the second largest real part determine the relaxation time r 
and the dynamical exponent z by the scaling behavior t ^ . Our finding is that z remains 
the same with the usual ASEP N = 2. The key to the result is that the spectrum of H in any 
multi-species particle sector includes that for an = 2 sector. We systematize such spectral 
inclusion relations by exploiting the poset structure of particle sectors. As a byproduct we 
find a duality, a global aspect in the spectrum of the Hamiltonian, which is new even at the 
Heisenberg point p = q. 

There arc some other multi-species models with different hopping rules such as ABC 
model [EKKM] and AHR model [AHR, RSS], etc. The recent paper [KN] says that the 
AHR model still has the dynamical exponent z = 3/2. The multi-species ASEP in this 
paper is a most standard generalization of the one-spccics ASEP allowing the application of 
the Bcthc ansatz. Let us comment on this point and related works to clarify the origin of 
integrability of the model. In [AR], a formulation by Hecke algebra was given for a wider 
class of stochastic models including reaction-diffusion systems. A Bethe ansatz treatment 
was presented in [AB]. As is well known [Ba], however, there underlies a two dimensional 
integrable vertex model behind a Bethe ansatz solvable Hamiltonian. In the present model, 
the relevant vertex model is a special case of the Perk-Schultz model (called "second class 
of models" in [PS]). It should be fair to say that the nested Bethe ansatz for the present 
model is originally due to Schultz [Sc]. For an account of the Perk-Schultz i?-matrix in the 
framework of multi-parameter quantum group, see [OY] and reference therein. 

The layout of the paper is as follows. In section 2, we introduce the multi-species ASEP 
Hamiltonian together with its basic properties. In section 3, we explain how the spectral 
gap responsible for the relaxation time is reduced to the one-species case and determine 
the dynamical exponent as in (3.17). Our argument is based on spectral inclusion property 
(3.11) and conjecture 3.1 supported by numerical analyses. In section 4, we elucidate a new 
duality of the spectrum of the Hamiltonian in theorem 4.12. An intriguing feature is that it 
emerges only by dealing with all the basic sectors of the N = L-state model on the length 
L ring. (The term "basic sector" will be defined in section 2.) The argument of section 4 is 
independent from the Bethe ansatz. In section 5, we discuss the Bethe ansatz integrability of 
the multi-species ASEP and the underlying vertex model including the completeness issue. 
We execute the nested algebraic Bethe ansatz in an arbitrary nesting order, which leads to 
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an alternative explanation of the spectral inclusion property. Except the arbitrariness of the 
nesting order, section 5.1.1 is a review. 

Appendix A gives a proof of the dimensional duality (theorem 4.9) needed in section 4. 
Appendix B is a sketch of a derivation of the stationary state by the Bethe ansatz. Appendix 
C contains the complete spectra of transfer matrix and Hamiltonian in the basic sectors 
with the corresponding Bethe roots for {p,q) = (§,5) and L = 4. They agree with the 
completeness conjectures in section 5.1.2. 



2 Multi-species ASEP 
2.1 Master equation 

Consider an L-site ring Zl where each site i G Z^, is assigned with a variable (local state) 
ki G {1, . . . , A^} {N > 1). We introduce a stochastic model on such that nearest neighbor 
pairs of local states (a,/3) = {ki, fci+i) are interchanged with the transition rate: 

fp iia>(3, 

ap^ Pa < (2.1) 
\q it a < p, 

where p and q are real nonnegative parameters. More precisely, the dynamics is formulated in 
terms of the continuous-time master equation on the probability of finding the configuration 
(fci, . . . , kz) at time t: 

^P{ki,...,kL;t) = ^ Q{ki+i - ki)P{ki,. . . ,ki-i,ki+i,ki,ki+2, ■ ■ ■ ,kL;t) 

i&^L ^2 2) 

- e{ki-ki+i)P{ki,...,kL;t), 



where 9 is a step function defined as 

(p (x > 0), 

G(x) = <j (x = 0), (2.3) 
[q (x<0). 

(Actually, 0(0) can be set to any value.) 

Our model can be regarded as an interacting multi-species particle system on the ring. 
We interpret the local state fcj = a as representing the site i occupied by a particle of the 
ath kind. The transition (2.1) is viewed as a local hopping process of particles. We identify 
the first kind particles with vacancies. Note that a particle a has been assumed to overtake 
any /3(< a) with the same rate p as the vacancy /3 = 1. 

We call this model the multi-species ASEP or more specifically the (A^— l)-species ASEP. 
The usual ASEP corresponds to N = 2. We will be formally concerned with the zero-species 
ASEP {N = 1) as well. The case p = q will be called the multi-species symmetric simple 
exclusion process (SSEP). 

Let |1 ),..., |A^) be the basis of the single-site space and represent a particle config- 
uration {ki, . . . , kL) as the ket vector . . . , ki) = \ki) (g) • • • (g) |A;l) G (C-^)®^. In terms of 



3 



the probability vector 



\P{t))= P{ki,...,kL;t)\ki,...,kL), (2.4) 

l<*;i<iV 

the master equation (2.2) is expressed as 

j^\P{t)) = H\P{t)), (2.5) 
where the hnear operator H has the form 

H=Y, hi+i, (2-6) 

h= ^ {-pEpii®Eaa-qEaa®Epi3+VEal3®Epa + qEiia®Eal}). (2.7) 
l<a<l3<N 

Here /ij^j+i acts on the ith and the (i + l)th components of the tensor product as h and as 
the identity elsewhere. Eap denotes the N hy N matrix unit sending I7) to (^-y^|a). 

The equation (2.5) has the form of the Schrodinger equation with imaginary time and 
thus provides our multi-species ASEP with a quantum Hamiltonian formalism. Of course in 
the present case, P{ki, . . . , /c/,; t) itself gives the probability distribution unlike the squared 
wave functions in the case of quantum mechanics. Nevertheless we call the matrix H the 
Hamiltonian} in this paper by the abuse of language. 



2.2 Basic properties of Hamiltonian 

Our Hamiltonian H is an by matrix whose off-diagonal elements are p, q or 0, and 
diagonal elements belong to pZ<o + qZ<Q. Each column of H sums up to assuring the 
conservation of the total probability X^i<fe.<jv Pi^i, ■ ■ ■ , ^l; t). It enjoys the symmetries 

[H, C] = 0, RHR-^ = QHQ-^ = H\p^q. (2.8) 

where C, R and Q are linear operators defined by 

C\ki,...,kL) = \kL,ki,...,kL-i) (cyclic shift), (2.9) 

R\ki,...,kL) = \kL,kL-i,...-,ki) (reflection), (2.10) 

Q\ki,...,kL) = \N+l-ki,...,N+l-kL) ("charge conjugation") (2.11) 

satisfying = R^ = = 1. H is Hermitian only at p = q, where it becomes the 
Hamiltonian of the ,s/(A^)-invariant Heisenberg spin chain H = P^i^ij^{(^i.i+i — 1) with (Tj^j+i 
being the transposition of the local states ki fci+i. In general, imaginary eigenvalues of H 
form complex conjugate pairs. In section 5, Bethe ansatz integrability of H for general p, q 
will be demonstrated. Although H is not normal in general, we expect that it is diagonalizable 
based on conjecture 5.1 and the remark following it. This fact will be used only in corollary 
4.8. 

more proper terminology is Markov matrix. 
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In view of the transition rule (2.1), our Hamiltonian H obviously preserves the number 
of particles of each kind. It follows that the space of states (C^)®^ splits into a direct sum 
of sectors and H has the block diagonal structure: 

(C^)®^ = 0F(m), V{m)= C\h,...,kL), (2.12) 

{fej} in sector m 

H = ^H{m), H{m) e'E^ndV{m). (2.13) 

m 

Here the direct sums 0^ extend over m = (mi, . . . , mjv) G such that mi + ■ ■ ■ + itin = 
L. The array m labels a sector V{m), which is a subspace of (C^)®^ specified by the 
multiplicities of particles of each kind. Namely, the latter sum in (2.12) is taken over all 
{ki, . . . , ^l) e {1, . . ■,N}^ such that 

mi m2 mjv 

SoTt{ki,...,kL) = 1^2^---'n^, (2.14) 

where Sort stands for the ordering non decreasing to the right. In other words, the {A;^} in 
(2.12) runs over all the permutations of the right hand side of (2.14). The array m itself will 
also be called a sector. A sector, m = (2, 0, 4, 1, 0, 3) for instance, will also be referred in terms 
of the Sort sequence as 1^3'^46^ = 1133334666. By the definition dimy(m) = ^I/Ili^i 

By now we have separated the master equation (2.5) into the ones in each sector ^l-P(i)) = 
H{m)\P{t)) with \P{t)) G V{m). Since the transition rule (2.1) only refers to the alternatives 
a > (5 or a < P, the master equation in the sector 1^3^46^ for example is equivalent to that 
in aP'h'^cd? for any l<a<h<c<d<N . 

Henceforth without loss of generality we take N = L and restrict our consideration to 
the basic sectors that have the form m = (mi, . . . , m„, 0, . . . , 0) for some n with mi, . . . , m„ 
being all positive. The basic sectors are labeled with the elements of the set 

M = {(mi,. . . ,m„) e Z>i I 1 < n < L, mi -I hm„ = L}. (2.15) 

In this convention, which will be employed in the rest of the paper except section 5, n plays the 
role of N in the sense that H{m) for m = (mi, . . . , m„) G 7V4 is equivalent to the Hamiltonian 
of the (n - l)-species ASEP. 

In section 3 we study specific eigenvalues of H that are relevant to the leading behavior 
of the relaxation. In section 4 wc elucidate a spectral duality, a new global aspect of the 
spectrum of H, which has escaped a notice in earlier works mostly devoted to the studies of 
the thermodynamic limit L — > oo under a fixed N. 

3 Relaxation to the stationary state 

3.1 General remarks 

The initial value problem of the master equation ^\P{t)) = H{m)\P(t)) with \P{t)) € V{m) 
in the sector m = (mi, . . . , m^) G M. (2.15) is formally solved as 

|P(t)) =e*-f^('")|P(0)). (3.1) 
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There is a unique stationary state corresponding to the zero eigenvalue of H{m). The sta- 
tionary state has a non-uniform probabiUty distribution [PEM] except for the zero-species, 
the one-species and the SSEP cases, where P{ki, . . . , k^) = 1/ dim V{m). All the other eigen- 
values of H{m) have strictly negative real parts, which are responsible for various relaxation 
modes to the stationary state. (The associated eigenvectors themselves are not physical prob- 
ability vectors having non-negative components.) Let us denote by |P(oo)) the stationary 
state. In general, the system exhibits the long time behavior 

|P(t)) - |P(oo)) ~ e-*/- (t^oo), (3.2) 

where r is the relaxation time. An important characteristic of the non-equilibrium dynamics 
is the scaling property of the relaxation time with respect to the system size: 

r ~ (L ^ oo), (3.3) 

where z is the dynamical exponent. The thermodynamic limit L ^ oo is to be taken under 
the fixed densities pj = rrij/L for j = 1, . . . , n. (Recall L = mi + ■ ■ ■ + rUn in (2.15), therefore 

Pi H hpn = !•) 

The spectrum in a sector m e M will be denoted by 

Spec(m) = multiset of eigenvalues ofH{m), (3-4) 

where the multiplicities counts the degrees of degeneracy. Spcc(m) is invariant under complex 
conjugation. The charge conjugation property (2.11) implies the symmetry 

Spec(mi, . . . , m„) = Spec(m„, . . . , mi). (3.5) 

We say that a complex eigenvalue x of H{m) is larger (smaller) than y if Re(x) > Rc(y) 
(Re(a;) < Re(y)). Spcc(m) contains as the unique largest eigenvalue. For a finite L, we say 
an eigenvalue E of H{m) is second largest if 

Re ^ = max Re(Spec(m) \ {0}) . (3.6) 

From (3.1) and (3.2), the scaling property (3.3) is equivalent to the following behavior of the 
second largest eigenvalues 

max Re(Spec(m) \{0}) = -cL-^ + o(L-^) (L ^ oo) (3.7) 

if the initial condition |-P(0)) is generic. Here c > is an "amplitude" which can depend on 
pi, . . . , Pn in general, but not on L. 

In the remainder of this section, we derive the exponent z of the multi-species ASEP 
based on (3.7). Our argument reduces the problem essentially to the one-species case and is 
partly based on a conjecture supported by numerical analyses. 

3.2 Known results on the one-species ASEP 

In this subsection we review the known results on the one-species ASEP. Thus wc shall 
exclusively consider the sector of the form m = (mi, 777,2) £ -M, and regard the local states 1 
and 2 as vacancies and the particles of one kind, respectively. Recall also that L = mi +m2- 
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Figure 1: Spec(3,4) for {p,q) = (0.8,0.2). 



The second largest eigenvalues are known to form a complex-conjugate pair, which will 
be denoted by E^{m). See figure 1. When mi = m2 or p = q, the degeneracy E~^{m) = 
£"{171) £ M occurs. 

In [GS, K, GM], the large L asymptotic form 

E^{{1 - p)L,pL) = ±27ri|(p - ^)(1 - 2p)\L-^ - 2C\p - qWp{l-p)L~l + 0{L-^) (3.8) 

with a fixed particle density p = m2/L was derived for p 7^ g by an analysis of the Bethe 
equation. The two terms are both invariant under p ^ ^ ~ P- The constant C has been 
numerically evaluated as C = 6.50918933794. . . . Thus from (3.7) the one-species ASEP for 
p q has the dynamical exponent 2; = |, which is a characteristic value for the Kardar- 
Parisi-Zhang universality class [KPZ]. 

In the SSEP case p = q, the Hamiltonian H{m\^m2) is Hermitian, hence all the eigen- 
values are real. The system relaxes to the equilibrium stationary state. For a finite L, the 
second largest eigenvalues take the simple form 

E'^ {mi,m2) = E~ (1711,1712) = — 4psin^(— ) (0 < 1712 < L), (3.9) 

which is independent of the density p = m2/L as long as < p < 1. The asymptotic behavior 
in L — > cxD is easily determined as 

^^((1 - p)L, pL) = -AiT^pL-'^ + 0(L-^), (3.10) 

3 

which is free from a contribution of order L~2 . From (3.10), we find the dynamical exponent 
z = 2, which is the characteristic value for the Edwards- Wilkinson universality class [EW]. 
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3.3 Eigenvalues of H{2, 1, 3, 1) -an example 



Let us proceed to the multi-species ASEP. Before considering a general sector in the next 
subsection, we illustrate characteristic features of the spectrum along an example. Figure 2 
(a) is a plot (black dots) of the spectrum Spcc(2, 1, 3, 1) on the complex plane. We recall 
that the sector (2, 1, 3, 1) means the ring of length 7 populated with 4 kinds of particles with 
multiplicities 2, 1,3 and 1, among which the first kind ones are regarded as vacancies. 

For comparison, we have also included the plot of the spectra in the one-species sec- 
tors (2,5), (3,4) and (6,1) in different colors and shapes. These one-spccics sectors are 
related to the multi-species sector (2, 1,3, 1) as follows. The sector (2,5) = (2, 1 + 3 -|- 1) 
is obtained by identification of all kinds of particles (except for vacancies) as one kind of 
particles. The sector (3,4) = (2 -|- 1, 3 -|- 1) is obtained by identification of the second 
kind particles as vacancies and the rest of particles as one kind of particles. The sector 
(6, 1) = (2 -|- 1 -h 3, 1) is obtained by identification of the second and the third kinds of parti- 
cles as vacancies. In figure 2 (a), we observe that all the colored dots overlap the black dots. 
Namely, Spec(2, 5), Spec(3, 4), Spec(6, 1) are totally embedded into Spec(2, 1, 3, 1). 

Figure 2 (b) shows that the second largest eigenvalues (denoted by E^{2, 1, 3, 1)) in those 
one-species sectors form a string within Spec(2, 1, 3, 1) near the origin. Although their real 
parts are not strictly the same, there is no black dot between the string and the origin. More 
precisely, there is no eigenvalue in the sector (2, 1, 3, 1) which is nonzero and larger than any 
second largest eigenvalues in the one-species sectors (2,5), (3,4) and (6, 1). This property is 
a key to our argument in the sequel. 

3.4 Eigenvalues of multi-species ASEP 

Let us systematize the observations made in the previous subsection. First we claim that the 
following inclusion relation holds generally: 

Spec(mi, . . . , rrin) D Spec(mi -|- • • • -|- mj,mj^i -|- • • • -|- m^) (1 < J < w — 1). (3.11) 

Each one-species sector appearing in the right-hand side of (3.11) is obtained by the identi- 
fication similar to the previous subsection: 



The relation (3.11) is a special case of the more general statement in theorem 4.5. See also 
section 5.2.1 for an account from the nested Bethe ansatz. 



Next we introduce a class of eigenvalues of H{m) for a multi-species sector m = (mi, . . . , m^) G 
M by 



Re(£^^(m)) = Re{EJ{m)), but the subscript j does not necessarily reflect the ordering of 
the eigenvalues with respect to their real parts. Generalizing the previous observation on 
figure 2, we make 




miH hnrij 



mj+iH hm, 
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Figure 2: (a) Spec(2, 1, 3, 1) (black dot), Spec(2,5) (red x), Spec(3,4) (blue square) and 
Spec(6, 1) (green triangle) with {p,q) = (0.8,0.2). (b) An enlarged view near the origin. 
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Conjecture 3.1. In any sector m = (mi, . . . , m„) G M, there is no eigenvalue E G Spec(m) 
such that 

max{Re£;j^(m),...,Re£;^_i(m)} <Re£;<0. (3.13) 

The one-species case n = 2 is trivially true by the definition. So far the conjecture has 

been checked in all the sectors m satisfying dimy(m) < 8000. 

Admitting the conjecture, we are able to claim that the second largest eigenvalues in 
Spec(m) are equal to E^{m) for some 1 < J < n — 1. (Such j may not be unique.) The 
asymptotic behavior of E^{m) is derived from (3.8) and (3.12) as 

Ef{piL,...,pnL) 

I 3 (3-14) 

= ±27ri|(p - g)(l - 2rj)\L-^ - 2C\p - q\^rj{l - rj)L-2 + 0{L-^) 

for p ^ q, where rj = pi + - ■ ■ + pj is fixed. We remark that the leading terms in (3.14) depend 
on j only through rj in the amplitudes. We call the eigenvalues Ef{m), . . . ,E^_^{m) next 
leading. Thus the second largest eigenvalues are next leading. All the next leading eigenvalues 
possess the same asymptotic behavior as the second largest ones up to the amplitudes as far 
as the first 2 leading terms in (3.14) are concerned. 

With regard to the SSEP case p = q, the stationary state is an equilibrium state. H{m) 
is Hermitian and E^{m) is real. We have the following explicit form as in the one-species 
case: 

Etim) = ■■■= E^_,im) = -Apsm\^). (3.15) 

In other words, the next leading eigenvalues E^{m) are degenerated in the SSEP limit p—q — > 
0. See figure 3.4, where the string of the next leading eigenvalues shrinks to a point on the 
real axis as p — g approaches 0. 

As the one-species case (3.10), we find 

Ef{p^L, . . . , p^L) = -47rV~' + 0{L-^) . (3.16) 

To summarize, the results (3.14) and (3.16) lead to the following behavior of the relaxation 
time r: 

ili for p^q 

r ~ <^ (L ^ oo). (3.17) 

I for p = q, 



Therefore we conclude that the dynamical exponent of the multi-species ASEP is independent 
of the number of species. It belongs to the KPZ universality class {z = |) for p ^ q, and to 
the EW universality class {z = 2) for p = q. 

We leave it as a future study to investigate the gap between the next leading eigenvalues 
and further smaller eigenvalues, which governs the pre-asymptotic behavior of the multi- 
species ASEP. 
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Figure 3: Degeneracy in the SSEP limit p — q ^ in the sector (2, 1, 3, 1). {p, q) is taken as 
(a) (0.8,0.2), (b) (0.7,0.3), (c) (0.6,0.4) and (d) (0.5,0.5). 



4 Duality in Spectrum 

Throughout this section, a sector means a basic sector as promised in section 2.2. We fix the 
number of sites in the ring L G Z>2. Our goal is to prove theorem 4.12, which exhibits a 
duality in the spectrum of Hamiltonian. 

4.1 Another label of sectors 

Set 

= {1,2,...,L-1}, (4.1) 
S = the power set of Q. (4-2) 

Recall that the sectors in the length L chain are labeled with the set Ai (2.15). We identify 
Ai with S by the one to one correspondence: 

M B m = (mi, . . . , m„) < — > {si < • • • < Sn-i} = s G 5 (4.3) 

specified via Sj = mi + m2 + • • • + rrij, namely, 

mi r?i2 m„ 

p"^'^'^ I o'"^''^''^ I ■ ■ ■ I o"^''^^'^. 



S2 



where the numbers of the symbols o and | are L and n — 1, respectively. For example the 
identification Ai <^ S for L = 4 is given as follows: 
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(1,1,1,1) ^{1,2,3} 
(1,1,2) ^{1,2} (1,2,1) ^{1,3} (2,1,1) ^{2,3} 

(1,3) ^{1} (2,2)<-.{2} (3,1) ^{3} 
(4) ^ 0. 

An element of 5 will also be called a sector. In the remainder of this section we will mostly 
work with the label S instead of M. There are tt<S = 2^~^ distinct sectors. We employ the 
notation: 

5 = \ s = complement sector of s. (4-4) 
For a sector s = {si < • • • < s„_i} G S, we introduce the set V{s) by (see (2.14)) 

-p^s) = {k = {ki,...,kL) I Sort(A;) = 2^ • • -'n-l.ln-l }, (4.5) 

where Sort stands for the ordering non decreasing to the right as in (2.14). For 5 = G <S, 
this definition should be understood as 7^(0) = {(1, . . . , 1)}. 
To each sector we associate the bra and ket vector spaces 

V:= C{k^,...,kLl ^.= C\k^,...,kL). (4.6) 

ker{s) ker{s) 

Here fci, . . . , fc^ G {1, . . . , L} stand for local states. For example if L = 3, one has 

V0 = C|111), 
V{i} = C|122) © C|212) © C|221), 
V{2} = C|112) © C|121) © C|211), 
Vn = V(i,2} = C|123) © C|132) © C|213) © C|231) © C|312) © C|321). 

Note that the vectors like |222) and |113) are not included in any Vg because we are concerned 
with basic sectors only. See (4.5). In general, one has 



L! 

Si!(s2-Si)! • • • (s„_i-s„_2)!(-^^-Sn-i)! 



dim Vg = dim Vg* = —^^ — — ^^7^- ^ (4.7) 



for s = {si < • • • < Sn-i} G S. 

Suppose 7W 9 m ■f-^> s G 5 under the correspondence (4.3). We renew the symbols H{m) 
and V{m) in (2.12)-(2.13) as 

Vg = V{m), Hg = H{m). (4.8) 

The set S is equipped with the natural poset (partially ordered set) structure with respect 
to C. The poset structure is encoded in the Hasse diagram [St], which is useful in our working 
below. In the present case, it is just the L — 1 dimensional hypercube, where each vertex 
corresponds to a sector. Sectors are so arranged that every edge of the hypercube becomes 
an arrow s — t meaning that s C t and [jt = [ts+ 1. There is the unique sink corresponding to 
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{1}^12 



3^11 



{1,2}^1233 



{l,2}f^^l23 



{l}-«-s-1222 



{1}-H>122 {2}<^>112 



0f^-lll ) 



{1,2,3}^1234 






{1,3)^1223 






{2}oll22 






0.«-s-llll 



{2,3}^1123 



{3}f^'1112 



(b) 








{2,4}-H-11223 










{4}^11112 



{l,4}-«-s-12223 



(d) 

Figure 4: Hasse diagrams (a), (b), (c) and (d) for L = 2,3, 4 and 5, respectively. Sectors are 
labeled by S (4.2) as well as the sequence Sort(A;) of local states as in (4.5). 



the maximal sector (7 G 5 and the unique source corresponding to the minimal sector S 5. 
See figure 4. 

We introduce the natural bilinear pairing between the bra and ket vectors by 



ei, • • • , ^lIji, • • • Jl) 



■6, 



kL,jL- 



(4.9) 



With respect to the pairing, V^* and Vi are dual if s = t and orthogonal if s 7^ t. 

Any linear operator G acting on ket vectors give rise to the unique linear operator G acting 
on bra vectors via I {ki, . . . , fc^lG I |ji, . . . , j^) = {ki, . . . , k^] I G\ji, . . . ,jL) I and vice versa. 



We write this quantity simply as {ki, . . . , . . . , j^) as usual, and omit and unless 

an emphasis is preferable. The transpose G'^ of G is defined by {k\G'^\j) = (j|G|A;) (=: Gjk) 

for any k = {ki, . . . jk^) and j = {ji, . . . , ji). Of course is equivalent to G in the sense 
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that 



(4.10) 



4.2 Operator (ps,t 

Let 5, t G 5 be sectors such that s C t. We introduce a C-Unear operator (pg^i in terms of its 
action on ket vectors (p^ ^ : Vi Vg. We define ipg^^ to be the identity operator for any sector 
5. Before giving the general definition of the case s C t, we illustrate it with the example 
s = {2, 5} C t = {2, 3, 5, 8} with L = 9. The Sort sequence of the local states in the sense of 
(4.5) for r{i) and r{s) read as follows: 



2 3 5 8 

V{t = {2, 3, 5, 8}) : 11 1 2\ 33 1 444| 5, 

2 5 

r{s = {2,5}): 11|222|3333. 



(4.11) 



According to these lists, we define fg^i to be the operator replacing the local states as 3 — »• 
2, 4 — 3, 5 — 3 (keeping 1 and 2 unchanged) within all the ket vectors |A;i, . . . , kL) in Vt. 

General definition of (ps,i is similar and goes as follows. Suppose t = {^i < • • • < tn} and 
6 = t\ {tjj ,t^}. Then (pg^i is a C-linear operator determined by its action on base vectors 
as follows: 



\ki,...,kL) 



1^1' • • • ; ^l)' 



(4.12) 



where x' = x — | ij < x}. 
Example 4.1. 



'Pl2,124\(l>) 
'PU,124\(I>) 



|21433) 
121333) 

121322) 



112343) eVt (t = {l,2,4}, L = 5), 
112333) ^12,12M = 121222) - |12222), 

112232), 



<^24,i24|<^) = 111322) - 111232), 



^1,12 9^12,1241 
^9^2,12 9^12,1241 

</5l,14 Vl4,124| 
^f4:,14 <^14,124| 

f2,2A V'24,124| 
<y?4,24 ¥^24,124 1 <;^, 



0, 



|21222) - |12222), 
|11211) - |11121), 



0, 

|11211) - |11121), 



9^1,1241 
9^4,1241 
^2,12i\ 



|21222) 
111211) 



|12222), 
|11121), 



= 0, 970,124 1 W = 0, 
where (^12,124 is an abbreviation of ¥'{i,2},{i,2,4}) ^tc. 



The following property of ips,t is a direct consequence of the definition. 
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Lemma 4.2. For a pair of sectors s d i, let Sq d Si d ■ ■ ■ d Si he any sectors such that 
5o = = t and (tSj+i = (1% + 1 for allO < j < I. Then, 

In particular, the composition in the right hand side is independent of the choice of the 
intermediate sectors Si, . . . ,S;_i. 

In example 4.1, one can observe, for instance, V'4,124|0) = 954.14 9514,12410) = ¥^4,24 ¥'24,124|0)- 
Let us turn to the transpose ipj^. By the definition (see (4.10)), we have 

\h,...,kL) ^ E'l-?!'---'-?'^)' ^^-^^^ 

i6P(t) 

where S' extends over those j = (ji, . . . , jl) € 7^(t) such that ips,t\ji, ■ ■ ■ , Jl) = l^i, ■ ■ ■ , ^l)- 
For example in example 4.1, one has 

¥'14,124 121322) = 121433) + |31423) + |31432). 
From (4.13) and (4.12) it follows that (/9g^t¥'^t — dimVs which actually means 

dimyt^, dimVt 
^^'^^^'^ =dimt^^^^- ^^'^^^'t =dh^rT^^^^^* 

for any sectors 5 C t. As a result, we obtain 

Lemma 4.3. Let s d i be any sectors. 

(1) (p^^:Vi^Vs is surjective. 

(2) ^g^t : V* — s- V^ is injective. 

The kernel of (pg^^ and the cokernel of ipg ^ will be the key in our derivation of the spectral 
duality in section 4.5. 

By now it should be clear that </7t\{n},t ^^^^ vectors in a sector t or send them to the 
neighboring smaller sector t \ {n} in the Hasse diagram against one of the arrows. Similarly, 
Vs,su{n} never kills bra vectors in a sector s and send them to the neighboring larger sector 
5 U {n} in the Hasse diagram along one of the arrows. 

4.3 Commutativity of cp^^i and Hamiltonian 

The action Hg : Vg Vg oi our Hamiltonian (4.8) is specified by (2.6)-(2.7) as 

Hg\ki,...,kL) = ^ Q{ki - ki+i){\ki...ki+i,ki...kL) - \ki . ..ki,ki+i . ..kL))- (4.15) 



Proposition 4.4. ipg^i is spectrum preserving. Namely, cpg^iHi = Hg (fg^i holds for any sectors 
Sdt 
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Proof. Consider the actions on the ket vector \k) = \ki, . . . ,kL) € Vt'. 



<p,,iH,\k) = - ki+i){\k[ . . . ki+„kl ... ^1) - \k[ ...k[, k',+, . . . k'^)) , (4.16) 

H,^,,,\k) = ^ eik', - k',+,){\k[ . . . k',+„kl ...k'j}- \k[ ...k[, k',+, . . . k'r^)), (4.17) 

where x' is the one specified in (4.12). For simpHcity, let us write {ki,kij^i) as {x,y). Prom 
(4.12), we see that x > y impUes x' > y', and similarly x <y implies x' < y' . Prom this fact 
and the definition of B in (2.3), the discrepancy of the coefficients Q{x — y) and B(x' — y') 
in the above two formulas can possibly make difference only when (x > y and x' = y') or 
{x < y and x' = y'). But in the both cases, the vector \...y' ,x' ...) — \...x' ,y' ...) is zero. Thus 
the right-hind sides of (4.16) and (4.17) are the same. □ 

Our Hamiltonian arises as an expansion coefficient of a commuting transfer matrix r(A) 
with respect to the spectral parameter A. See (5.8). However, the commutativity ips,iT{\)i = 
r(A)gV3s^t does not hold in general. 

To each sector 5 = {s\ < ■ ■ ■ < s„_i} G S, we associate 

Spec(s) = multiset of eigenvalues of Hg, (4-18) 

where the multiplicity of an element represents, of course, the degree of its degeneracy. This 
definition is just a translation of (3.4) into the notation (4.8). The property (3.5) reads 

Spec(si, . . . , Sn-i) = Spec(L - . . . , L - si). (4.19) 

One has (JSpec (5) = dim 14 = dimV*. Lemma 4.3 (2) and proposition 4.4 lead to 

Theorem 4.5. There is an embedding of the spectrum Spec(5) ^ Spec(t) for any pair of 
sectors such that s C i. In particular, Spec(f2) contains the eigenvalues of the Hamiltonian 
Hg of all the sectors s G 5. 

See figure 5 for example. 



4.4 Spectral duality in the maximal sector Q, 

As indicated in theorem 4.5, the structure of the spectrum in the maximal sector O G 5 is of 
basic importance. In this subsection we concentrate on this sector and elucidate a duality. 
Define a C-linear map lo by 

^- ' y» u 20) 

{ki,...,kL\ ^ sgn{k) \kL,. . . ,ki), 

where sgn(fc) = sgn(A;i, . . . , A;l) stands for the signature of the permutation. (Note that V{^) 
is the set of permutations of (1,2,..., L).) Obviously, oj is bijective. 

It turns out that uj interchanges the eigenvalues of Hamiltonian as E —L{p + q) — E. 

Theorem 4.6. Let {(j)\ G be an eigenvector such that {4>\Hfi = E{(l)\. Set \^) = oo{{4>\) G 
Vn. Then Hn\il^) = {-L{p + q) - E)\tP) holds. 
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^ {l,2,3}ol234 

-4u, -^u-iv, -Su-iv, -3^, -3u, -3u, -3if + iiJ, 



- 3u + ii;, - 



, -"I + ^ r, -2ii, -2?f, -2?/, 
-u-\-iv, —u-\-iv, 



{l,2}^1233 

-3u-ii', -3u, -3ii + i7;, 



I 



-2^, -2u, -4 u- 



r, -7(-ri;, 
u + iv, 



{l,3}ol223 

-3ii, -3^, --J u-^ r, 



~T ^ + T ^' -2u, -2u, 







{l}^1222 ^ 

-2u, -u-iv, -u-\-iv, 



-u-iv, -u-iv, -u, 
-u-\-iv, -'U + ri;, 



{2,3}oll23 

-Su-iv, -Zu, -Zu-\-iv, 



-u, -u-\-iv. 



{2}-H-1122 



- z( + 4- r, -u, 







{3}f^lll2 

-2iL, —u-iv, -u-\-iv, 



0-H-llll 





Figure 5: Spec(s) for L = 4. u = p + q, v 
(4.19) can be also observed. 



q, r = \J —Ip^ + ISpg — 7(7^. The symmetry 



Proo/. Let ((/)| = XlfcG-PCn) /(^i' • • • , • • • , Then ((/"l-ffn = is expressed as 

Qih - k+i){fik^'^) - f{k)) = Efik), 

where we have used the shorthand k = {ki, . . . ,ki, fcj+i, . . . , k^) and feW = (ki, . . . , /cj+i, ki, . . . , 
Adding (p + q)Lf{k) to the both sides we get 

^ eih - h+^)f{k^^y) +Y^{p + q- eih - h+i))f{k) = {E + L{p + q))fik), 

Since kiS are all distinct in the sector VI under consideration, the coefficient in the second 
term equals 0(A:j+i — ki). Multiplication of — sgn(A:) = sgn(A;(*)) on the both sides leads to 

&iki - A:,+i)sgn(A:«)/(fe«) - ^ ©(^^+1 " ki)sgnik) f (k) 

= {-E-L{p + q))sgn{k)f{k). 

This coincides with the equation Hi^\ip) = {—L{p + q) — E)\Tp) on 

= Ukarm sgn(A;)/(/ci, . . . , kL)\kL, ...,ki). □ 

Remark 4.7. It is easy to see that {(p\ = J2keV(n)(ki, • • • > S is the eigen bra vector 
with the largest eigenvalue E = 0. It follows that uj{{(j)\) G is the eigen ket vector with 
the smallest eigenvalue —L{p + q). Namely, one has 

{H^ + L{p + q)) Y sgn{k)\kL,...,ki) = 0. (4.21) 
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In view of conjecture 5.1 and the remark following it, we assume the diagonalizability of 
the Hamiltonian Hq,^ . Then every eigenvalue in Spec(r2) is associated with an eigenvector 
in V^. Therefore theorem 4.6 implies 



Corollary 4.8. Spec(17) 



L{p + q) — Spec(n) . 



Im 



2 



• •• • • ••• 



1 




,•* X,^, 



-6 



Re 



• •• • • ■•• 



- 1 



-2 



Figure 6: Spec(r2) for L = 6 and (p, = (0.8,0.2). The symmetry with respect to —L{p + 
q)/2 = —3 can be observed. 

Figure 6 is a plot showing this property. The property of interchanging the eigenvalues 
of Hamiltonian E —L{p + q) — E will be referred as spectrum reversing. Our main task 
in the sequel is to extend u; to a spectrum reversing operator between general sectors, and 
to identify the "genuine components" that are in bijective correspondence thereunder. This 
will be achieved as uj° in theorem 4.12. 

4.5 Genuine components X* and Yg 

Theorem 4.5 motivates us to classify the eigenvalues Spec (s) in a sector s into two kinds. One 
is those coming from the smaller sectors u C s through the embedding Spec(u) ^ Spec(5). 
The other is the genuine eigenvalues that are born at 5 without such an origin. Having this 
feature in mind we introduce a quotient X* of V* and a subspace of as 



We call X* and Yj the genuine component of V* and V^, respectively. (We set X^ = Vg* = 
C(l, . . . , 1| and l0 = V0 = C|l, . . . , 1).) The Hamiltonian Hg acts on each X* and Yg owing 

^Theorem 4.5 is derived on the basis of generalized eigenvectors hence its vahdity is independent of the 
diagonalizabihty of the Hamiltonian. 




(4.22) 



18 



to proposition 4.4. The vector spaces X* and Yg are dual to each other canonically, therefore 

dimX* =dimy3. (4.23) 

We wish to focus on the spectrum that are left after excluding the embedding structure 
explained above and in theorem 4.5. This leads us to define the set of genuine eigenvalues of 
a sector s as 

Spec°(s) = multiset of eigenvalues of Hg\x* 
= multiset of eigenvalues of Hg\Y^. 

Let us write the image of (01 G V* in X* under the natural projection by [(0|]. Fix an 

embedding of X* into sending each eigenvector [{(f>\] G X* to an eigenvector G with 
the same eigenvalue satisfying [{<i>\] = [{<P'\]- The image of the embedding is complementary 
to X^m^^Im iPu,s^ therefore we can treat the first relation in (4.22) as V* = X* (BJ2u(zs^^ ^u,s- 
Then the following decomposition holds: 

^; = ©^:k,.- (4.25) 

uCs 

Prom theorem 4.5 and (4.25) we have 

Spec(s) = IJ Spec°(u), (4.26) 

uCfi 

where the multiplicity is taken into account for the union of the multisets. In terms of the 
cardinality, this amounts to 

dim Vg* = dim X* . (4.27) 

uCs 

Theorem 4.9 (Dimensional duality). For any sector s G S, the following equality is valid: 

dimX* = dimX|, 

or equivalently ttSpec°(5) = ttSpcc°(5). Here s denotes the complement sector (4.4). 

See figure 7 for example with L = 4. The proof is due to the standard Mobius inversion 
in the poset S and available in appendix A. 

The following lemma, although slightly technical, plays a key role in our subsequent 
argument. 

Lemma 4.10. 

(V ^n\{r},n'^(Im^n\{r},n) = for any reft. 
(^) ^s,n t^(Im^u,n) = unless u D s. 

Proof. (1) For brevity we write = \ {r}. We illustrate an example L = 5, = 
{1, 2, 3, 4}, = {1, 3, 4}, from which the general case is easily understood. Recall the scheme 
as in (4.11): 

12 3 4 

7'(J^ = {1,2,3,4}) : 1|2|3|4|5, 
P(J^2 = {1,3,4}) : lj2 2|3|4. 
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'{l,2,3}ol234 
I (24,1) I 




Figure 7: The data (dim 1^*, dim X*) = (tjSpec(s), ttSpec°(s)) is presented for each s in the 
same Hasse diagram (c) in figure 4. The dimensional duahty (theorem 4.9) can be observed. 
For a systematic calculation of these data, see appendix A. 



Thus n is the operator replacing the local states 3 — > 4, 4 ^ 5 and moreover changes 
(...2, 2, ...| into the symmetric sum (...3, 2, ...| + (...2, 3, ...|. At the next stage, to in 
(4.20) attaches the factor sgn(A;) which makes the above sum antisymmetric. Finally, ^p^^ ^ 
makes the antisymmetrized letters 2 and 3 merge into 2 again (and also does 4 — > 3, 5 — > 4), 
which therefore kills the vector. For example, 

(423121 'C^" (524131 + (53412| 
^ -131425) + 121435) 

^ -|21324) + 121324) = 0. 

(2) Note that Imc/j^f, = V*Lpy^,^. Thus we are to ask when fs,n^{^u^u,a) vanishes. It is 
helpful to view this as a process in the Hasse diagram going from V* to Vg via the maximal 
sector as in figure 8, where u = {ui, . . . , and s = {si, . . . ,5;,}. In figure 8, the arrows 
represent the factorization ^p^^q^ = ¥'u,uu{«a} " " " due to lemma 4.2 growing u up to 

Q. by adding Hj's one by one. Similarly the arrows \ stand for ip^^^ = v^g^jsulsb} ' ' ' ¥'n\{si},r! 
shrinking 17 down to 5 by removing Sj's one by one. (The arrows attached to Ui (si) are the 
same (opposite) as those in the Hasse diagram.) In this way 

^a,fit^(KVu,f2) = • • • ^n\{si},n^{- ' ' ^n\{ui},(i) 

= ■ ■ ■^n\{s,},n^{- ■ ■'^n\{u,},n) for any 1 < i < 6, l<j<a, 

where the second equality is due to lemma 4.2 which assures that the factorization is possible 
in arbitrary orders. From the assertion (1) we thus find that this vanishes if 5 Pi u 7^ 0. In 
other words, ip^ ll!{V* ip^^ n) = unless = 5 H u, or equivalently u 5 5. □ 
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a; 



ui 



5 



Figure 8: A conceptual scheme of the proof of lemma 4.10 (2). 



Proposition 4.11. 

Proof. (1) 



LemA3(l) 
— 



V^^ /'y*^ \ Lem.4.10(2) ^ ^ / 



n 
u3s 



(4.28) 



Taking the dimensions, we have 

dim Vs = dim ^ „ cj (x* ^„ j < ^ dim w (x* ^„ ^ 



Th.4.9 T ^* (4-27) 



< ^dimX* = ^ J^dimX^ = ^dimX^ = ' dimF/ = dim^. 

Thus all the inequalities < here are actually the equality =. Moreover, all the sums Yl 
(4.28) must be the direct sum ©, finishing the proof. 

(2) Let Yg = 99 J ^{^^ n)- argument similar to the proof of lemma 4.10 (2), one 

can easily show that Yg is killed by Vs\{n},s £ s. In view of lemma 4.2, this implies 

^ The proof is finished by noting dimYg = dim ip^ ^u>{X^ ip^ ^) = dimX| "^^-^ 
dimX; ^^=^' dimYg. □ 

Combining proposition 4.11 (2) and theorem 4.6, we arrive at our main result in this 
section. 

Theorem 4.12 (Spectral duality). For any sector s E 5 and its complementary sector s, 
there is a spectrum reversing bijection iO° between their genuine components: 

^ (4.29) 
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In particular, the genuine spectrum enjoys the following duality: 



Spec°(5) = -L{p + q)- Spec°(s). 



(4.30) 



This relation is a refinement of theorem 4.9. 

Example 4.13. Figm'e 9 presents Spec°(5) for L = 4 in the same format as figure 5. All the 
genuine eigenvalues form pairs with those in the complementary sectors to add up to —L{p + 
q) = —Au including the multiplicity. The full spectrum Spec(s) in figure 5 is reproduced from 
the data in figure 9 and (4.26). 



{1,2,3}^1234 



{1,2}^1233 

-Zu-iv, -Zu-\-w, -2u 




{1}^1222 

-2u, -u-iv, -u-\-iv 



{1,3)^1223 

-y r, -u 

{2}^1122 



-3u, -4 u-4r r, 



{2,3}^1123 

- Su - iv, - Su + itJ, - 2u 




{3}^1112 

2u, -u-iv, -u-\-iv 



0^1111 





Figure 9: Spec°(s) for L = 4. u = p + q, v = p — q, r = ■\/ —Ip^ + \%pq — 7g^. 



Remark 4.14. The genuine spectrum Spec° also enjoys the symmetry (4.19). It follows that 
if a sector t satisfies t D s,s with s = (si, . . . ,s„_i) and s = (L — ... ,L — si), then Hi 
is degenerated because of Spec(t) D Spec°(s) U Spec°(s) and Spec°(s) = Spec°(s). 

5 Integrability of the model 

Our multi-species ASEP is integrable in the sense that the eigenvalue formula of the Hamil- 
tonian can be derived by a nested Bethe ansatz [Sc]. See also [AB, BDV]. 

As mentioned in section 1, our Hamiltonian is associated with the transfer matrix of the 
Perk-Schultz vertex model [PS]. In section 5.1, we derive the eigenvalues of the transfer 
matrix in a slightly more general way than [Sc]. Namely we execute the nested Bethe ansatz 
in an arbitrary "nesting order". In section 5.2, we utilize it to give an alternative account of 
the spectral inclusion property (theorem 4.5) in the Bethe ansatz framework. We also recall 
the original derivation of the asymptotic form of the spectrum following [K]. In section 5.3, 
the Bethe ansatz results are presented in a more conventional parameterization with the 
spectral parameter having a difference property. 
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5.1 Nested algebraic Bethe ansatz 

5.1.1 Transfer matrix and eigenvalue formula 

Let us derive the eigenvalues of the Hamiltonian H (2.6) for the (A'^ — l)-species ASEP on 
the ring Z^, by using the nested algebraic Bethe ansatz. Let Wj be a vector space W = 
at the jth site of the ring. We define a matrix RjkW € End(VFj (g) Wk) as 

RjkW = Pjk{i + Xhjk), (5.1) 

where Pj^ and hj^ are, respectively, the permutation operator and the local Hamiltonian 
(2.7) acting non-trivially on Wj W^- The non-zero elements are explicitly given by 

RZW = 1, ^«2(A) = 1'^ " ^ ^' i?S(A) = 1'-'^ " < ^' (5.2) 
''"^ ^ '"^^ ^ \pX fov a>l3, '"^^ ^ \l-pA fora>/3. ^ ^ 

Here a,/3 e {1, 2, . . . , iV}, and Rl^^iX) stands for R,k{X){\a)j ^ |/3)fe) = h)j ® |5)fc^5(A) 
(summation over repeated indices will always be assumed). The above i?-matrix satisfies the 
Yang-Baxter equation [Ba] 

i?23(A2)i?13(Al)i?12(A) = i?12(A)i?13(Al)i?23(A2), (5.3) 

where the parameter A is given by 

A = ^(Ai, A2) = 1 , AT^I X X • (5-4) 
1 - (p + g)A2 +pgAiA2 

This is not a simple difference Ai — A2. However, one can restore the difference property by 
changing variables as in section 5.3. Thanks to (5.3), the transfer matrix T(A) G End(M^®''") 

T{\) = tTwo[RoL{\)---Roi{\)] (5.5) 
constitutes a one-parameter commuting family 

[r(Ai),r(A2)] = o. (5.6) 

It means that T{X) is a generating function for a set of mutually commuting "quantum 
integrals of motion" Ij {j = 1,2,...): 



I. = (|)'l..T(A) 



(5.7) 

A=0 



Iq is the momentum operator related to the shift operator C (2.9) by C = expXo- ^1 yields 
the ASEP Hamiltonian H (2.6): 



Thus the eigenvalue problem of H is contained in that of T(A). To find the eigenvalues of 
r(A), we introduce the monodromy matrix T(A) G End(VFo ® W®^) by 

r(A) = i?oL(A)---i?oi(A). (5.9) 
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Its trace over the auxiliary space Wq reproduces the transfer matrix (5.5) 

r(A) = tTwoT{X). 
Prom the Yang-Baxter equation (5.3), one sees the following is valid: 

r2(A2)ri(Ai)i?i2(A) = i?i2(A)ri(Ai)T2(A2), 



(5.10) 



(5.11) 



where i?i2(A) here acts on the tensor product of two auxiliary spaces. 

Let us define the elements of the monodromy matrix in the auxiliary space as T(A)|a)o 
Ta{X)\P)o, where t£{X) acts on the quantum space W^^. More explicitly, 



r(A) 



fT-HX) Ba,{X) 
C-'^iX) T-^X) 



BaAX)\ 
-4"^ (A) 



\c--{x) r„«-(A) ••• r„»-(A)y 

Ba, (A) := (A), C"^- (A) := Tj^^ (A) for je{2,...,N}, 



(5.12) 



Here we have introduced the indices ai,...,aN that are arbitrary as long as {ojjjLi = 
{1, . . . , N}. They specify the nesting order |ai), {02), ■ ■ ■ , \aN) ^. 

Let |vac) := |ai)i (g) |ai)2 • • • <8) \ai)L be the "vacuum state" in the quantum space. It 
immediately follows that the action of T{X) on |vac) is given by 



r (A) I vac) 



/I Ba,{X) 

d(A)(g/p)^^i2 
\0 



^a.(A) \ 



diX){q/p)^0^^ ) 



I vac). 



(5.13) 



where 




for Oj > aj. 



d{X) := (pA)^, Oij := 9{ai - a,) : 
Using the relation (5.11), we can verify the following commutation relations: 

B^{X)Bp{X') - 



(5.14) 



r:i{x)B^{)^) 

T^^iX)B^iX') -- 



B0{X')B^{X) fora = /?, 

^5(e(A, X'))Bs{X')B^{X) for a^(5, 

f/(A',A)5,(A')r„»i(A) +5(A, A')S«(A)r„«i(A') for ai < «, 
[7(A', A)5„(A')r„«nA) +5(A,A')S«(A)r„«i(A') for ai > a, 

/(A, AO<JC(A, X'))Bs{>^)Ti{X) - g{X, X')B^{X)T£{X') for ai < (3, 



7(A, X')mi{aX, X'))Bs{X')Tf{X) - g{X, X')B^{X)Ta'{X') for ai > /3. 



(5.15) 



'in the standard nested algebraic Bethe ansatz, the nesting order is chosen as aj = j. 
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Here a,P,j,S G {cij}j=2' functions /, g, f and g are defined by 



/(A,/x) 



1 - {p + q)n + pqX/j, 



5(A,m) = 1-/(/x,A) = /(A,//)- 



]9(A - fl) 

q _ {l-qX)il-pn) 



= /(A,/x)|p^q = -/(A,/x), 5(A,/x) =£?(A,/x)|p^g. (5.16) 
Consider the following state with the number of particles of the a^-th kind being niaj 
|{AW}) = F°i-"iS,,(AS'^)---S„„^(A«)|vac), e {a2, . . . , oat} (1 < i < m), (5.17) 



where 



N 

rik ■■= ^ {1 < k < N - 1; no = L). 

j=k+l 



(5.18) 



The sum over repeated indices in (5.17) is restricted by the condition 

<j< ni, aj = Gk} = ma, {2 < k < N). (5.19) 
Then the action of T{X) on |{A(^)}) is calculated by using the relations (5.15) and (5.13): 

ni ni 



r(A)|{AW}) = 



N 



a=2 



ni 



+ i^"^---"-r('^£::::±\(A|{A«})d(A) n/(A, A«) (5-20) 

j=l j=l 

where the product YYjLi ^Pji'^f^) is ordered from left to right with increasing j; aj,Pj G 
{02, . . . , ajv}; rik is an integer given by 



N 

rik-- rUa^Okj il<k<N-l), 
j=k+l 



(5.21) 



and T^^^i\';::fal\ (A|{A(i)}) is a matrix element of T(i)(A|{A(i)}) G End((C^-i)®^) defined by 

(L-ni)6ia 



N 



tW(A|{aW}) = ^ 



a=2 



P 



^oni(^(A,A«))---i?oi(?(A,Ai'))) 



(5.22) 



Note that the sum corresponds to the trace over the {N — l)-dimensional auxiliary space 
spanned by the basis vectors |aj)o (2 < j < A^), and the quantum space acted on by 
r(^)(A|{A(^)}) is spanned by the vector (8)^j^^|aj)j where aj G {ak}k=2- 
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If we set F'^i" °="i as the elements of the eigenstate for T'^^)(A|{A(^^}) and choose the 
set of unknown numbers {A^-^^}J1^ so that the unwanted terms (u.t.) in (5.20) become zero 
(u.t. = 0), the eigenvalue, written A(A), of the transfer matrix T(A) is expressed as 

(\ ni "1 rii 
- n /(^r- ^) + A(^)(A|{A(i)})d(A) n /(A, Af). (5.23) 
^/ j=i j=i 

Here A(-'^)(A|{A'^^)}) is the eigenvalue of r(^)(A|{A'^^)}), which will be determined below. Noting 
that 

C(e(Ai,/i),e(A2,M))=aAi,A2), (5.24) 

and using the Yang-Baxter equation (5.3), one finds that the transfer matrix r(-'^)(A|{A'^^)}) 
forms a commuting family 

[r(Ai|{AW}),T(A2|{AW})]=0. (5.25) 

Hence the method similar to the above is also applicable to the eigenvalue problem of 
r(^^(A|{A(^)}). Namely constructing the state 



n2 



|{A(2)}) =F««i--"24i)(Af))...i?« (A(,2))|vac«), |vac«) ■.= (^\a,)j, (5.26) 

where aj G {3, . . . , N} and 

5W(A) = kni(e(A,A«))---i?oi(^(A,A«))l"' , (5.27) 



we obtain 



/ \ (L—ni)6i2 / \ n2 "2 

A<"(AKA»'})-(^) (^) n/(Af.A) 



3=1 

ni ^ n2 



+ A(^)(A|{A(^)})n— \(iy-nMAf). (5.28) 
j=i J l'^, '^j ) j=i 

Note that the coefficient i^'(i)°=i"'""2 in (5.26) and A*^^) (A|{A(^)}) in the above are, respectively, 
the elements of the eigenstate and the eigenvalue for the transfer matrix 

T(2)(A|{A(2)}) = ^ r-) kn2(€(A,A(^)))...i?.,(4(A,Af))) " . 



(5.29) 

The sum corresponds to the trace over the [N — 2)-dimensional auxiliary space spanned by 
the basis vectors |aj)o (3 < j < A^), and the space acted on by r(^)(A|{A(^)}) is spanned by 
^''j=iWj)j where aj G {ofcl^La- Repeating this procedure, one obtains 

A(0(Ai{A(0})=(^j n/(^r^^) 

n-i ^ rii+i 

+ A('+^)(A|{A('+i)})n-^ n /(A^A^)), (5.30) 
j=i /(A, Xj ) j=i 
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for 2 < / < iV - 2, and 
aM(A|{A(^-2)}) 



njv-1 



E,=i («j-i-%)fjiv-i-njv-r , 

n /(Af - 



+ 



p 



"JV-l 



n 



j=l /(A, A^- ^; 



n 



(5.31) 



Thus we finally arrive at the eigenvalue formula of the transfer matrix: 

-ni ni 



A(A) 



P 



^2 / X Ei=i(nj-i-nj)ejfe+i-nfe+i "fc 



fe=i 



P 



"fc+1 



n /(^'^. 



(iV-l)x 

'j > 



(5.32) 



The unwanted terms disappear when the set of unknown numbers {A|"^}"i^-|^ (1 < ^ < -/V — 
1) satisfy the following Bethe equations, which are also derived by imposing the pole free 
conditions on the eigenvalue formula: 



P 



(L-ni)6li2 



1 



fc=i /(Aj ,Aj^. )k=if{\ 



g\ Ej=iK-i-%)^ji+i-Ej=\K-i-%)^ji 



-^^+^^+1 /(Ai'\Af)nfcL7/(Af,Ag-^^) 

M/(Af,A«)nL"iV(Aruf) 



(2<Z<iV-2), 



\ E^i^ ("j - 1 )^'jJV -EjLi^ {nj-i-nj)ejN-i 



pj 



p 



.(AT-l) ,(JV-1)^ 11 •'V^i ''^fe 



fc=l /(.A- ,Aj. j 



(5.33) 



Inserting the expression (5.32) into (5.8), one finds the spectrum of the Hamiltonian: 



£ = ^lnA(A) 



(l-p\f){l-q\f) 



E 



A=0 ,=1 



A 



(1) 



(5.34) 



Though the explicit form of the Bethe equation (5.33) and its solutions depend on the 
nesting order in general, the spectrum of the Hamiltonian (5.34), of course, does not depend 
on it. 
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5.1.2 Completeness of the Bethe ansatz 

In this sub-subsection we exclusively consider the standard nesting order aj = j {1 < j < N). 
Let us recall our setting and definitions. We consider the transfer matrix T(A) (5.5) acting 
on the sector V{m). See (2.12). The data m = {mi, . . . jm^r) G specifics the number 
nij of the particles of the jth kind and rrii + ■ ■ ■ + rriN = L. The sector m is basic if it has 
the form m = (mi, . . . , nin, 0, . . . , 0) with mi, . . . , nin all positive for some 1 < n < L. The 
basic sectors are labeled either with M (2.15) or S (4.2) by the one to one correspondence 
M 3 m ^ 5 £ S (4.3). For a basic sector m, we write V{m) also as Vg as in (4.8). Yg is the 
genuine component of Vg defined in (4.22). Spec°(s) is the multiset of genuine eigenvalues of 
H in Ys (4.24). 

Conjecture 5.1. Suppose that p ^ q are generic. Then for any sector V(m) which is not 
necessarily basic, there exist d = dimF(m) distinct polynomials Ai(A), . . . ,K^(\) in A such 
that det(C-T(A)) = n'=i(C-A,(A)). 

We call Ai(A), A2(A), . . . , Ad(A) the eig en-polynomials of r(A). (It should not be confused 
with the characteristic polynomial dct(C — T(A)).) A direct consequence of conjecture 5.1 
is that the transfer matrix r(A) hence the Hamiltonian H are diagonalizable in arbitrary 
sectors. (At p = q, the diagonalizability still holds but Ag(A)'s are no longer distinct due to 
degeneracy caused by sZ(A?^)-invariance.) 

Now we turn to the completeness of the Bethe ansatz. In the remainder of this subsection 
and section 5.2.1, by the Bethe equations we mean the polynomial equations on {A^''' | 1 < 
I < N — 1,1 < j < Uj} obtained from (5.33) by multiplying a polynomial in them so that the 
resulting two sides do not share a nontrivial common factor. We say that a set of complex 
numbers {A^'^} is a Bethe root if it satisfies the Bethe equations. Bethe roots {A^'^} and 

{A^'^'} are identified if A^'^ = X^^^^' for some permutation ki,...,kni of 1, . . . , n; for each I. We 

say that a Bethe root {A^'^} is regular if none of them is equal to 1/p and two sides of any 
Bethe equation are nonzero. Using the same notation as in conjecture 5.1, we propose 

Conjecture 5.2 (Completeness). Suppose pi^ q are generic. 

(1) For any sector, all the eigen-polynomials Ap(A) are expressed in the form (5.32) in 
terms of some Bethe root. 

(2) For a basic sectors G S, there exist exactly dim Yj regular Bethe roots and the associated 
dirnY^ eigen-polynomials among Ai(A), . . . , Ac;(A). 

(3) The dirnY^ eigen-polynomials in (2) give Spec°(0) hy the logarithmic derivative (5.34). 

In view of section 4.5, it is natural to call the (conjectural) dim Yg eigen-polynomials in 

conjecture 5.2 (2) the genuine eigen-polynomials of the basic sector s. Then conjecture 5.2 
(3) is rephrased as claiming that the spectrum Spec(s) and the genuine spectrum Spec°(s) 
are obtained by the logarithmic derivatives of the dim^ eigen-polynomials and the diml^ 
genuine eigen-polynomials, respectively. 

Some examples supporting conjectures 5.1 and 5.2 are presented in appendix C. In 
conjecture 5.2 (1), the Bethe roots corresponding to a non-genuine eigen-polynomial are not 
necessarily unique. See the 2nd and the 3rd examples from the last in appendix C. We expect 
that the Bethe vectors associated with the regular Bethe roots form a basis of Y^. 
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Theorem 4.5 and conjecture 5.2 (2) bear some analogy with the sZ(iV)-invariant Heisenberg 
chain (p = q). There, the number of the Bethe roots are conjecturahy the Kostka numbers 
[KKR] and the spectral embedding is induced by sl{N) actions. Here, the analogous roles 
are played by dim Yg and t, respectively. 

5.2 Properties of the spectrum 

Now we derive some consequences of the eigenvalue formula (5.32), (5.33) and (5.34). Sections 
5.2.3 and 5.2.4 are reviews of known derivation for reader's convenience. 



5.2.1 Spectral inclusion property 

First we rederive the spectral inclusion property (theorem 4.5) in the Bethe ansatz framework. 
Consider the sector t G 5 where the number of particles of the j'th kind is rrij > 1 for any j. 
In the notation (4.3), t reads 

t = {mi, mi + m2, . . . ,rni + m2 H h mN-i} 

^ 1"'! 2™2 . . . iV^iv = . . . Q^-iv ^ (^5 35^ 

where mi + • • • + mjv = L>2. Set 



ajv = aN-i + 1, ^ 7 (1 < j < riN-i)- (5.36) 

Due to the relations 



1 



J p 



UN-i = 0, f[^,\)=Pf (x, 1") = 1 for A 7^ -, (5.37) 



\P J Q. \ P J P 
the following reduction relation holds: 

_ / ^.\'Efji(.iij-i-'n'jWjN+nN-i 

A(A) = A(A) + d{\) i^j . (5.38) 

Here A(A) stands for the eigenvalue formula of the (N — 2)-species case in the sector 

"l ■■■"jv— 1 — ' ' ' "iV— 1 "at ■ ■ ■ 1) 

^t\{miH l-majv_J = u, (5.39) 

where 

i, = h (l<j<N-l). (5.40) 

[Uj — 1 for ttj > Oat 

If {x'-p \ I < j < ni, 1 < I < N - 1} is a solution of the Bethe equation in the sector t, 

so is {A^'^ I l<j<ni, 1<1<N — 2} left after the substitution (5.36) in the sector u. 
This is because the last Bethe equation in (5.33) for (N — l)-species case becomes trivial, or 
alternatively one may say that the resulting (N — 2)-species Bethe equation guarantees that 
A(A) is pole-free. 
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Inserting (5.38) into (5.34), and using d{0) = d'{0) = 0, one thus finds the set of eigenval- 
ues of the Hamiltonian for the sector t includes that for the sector u. Applying this argument 
repeatedly, one can see 

Spec(s) C Spec(t) for set. (5.41) 

That is theorem 4.5. Since the solutions of the Bethe equation (5.33) depend on the nesting 
order, the set of solutions characterizing the above Spec(s) are, in general, not included in 
the original set of solutions characterizing Spec(t). 



5.2.2 Stationary state 

One of the direct consequences of the above property is that the stationary state E = for 
an arbitrary sector t (5.35) is given by setting all the Bethe roots to 1/p, i.e. 

A?^ ^ - {l<j<ni,l<l<N). (5.42) 

It immediately follows that the eigen-polynomial of the stationary state is given by 

A(A) = 1 + ^^ d(A), (5.43) 

k=l ^^'^ 

where is defined by (5.18), and we consider the standard nesting order aj = j {1 < j < N). 

On the other hand, in the framework of the Bethe ansatz, the calculation of the corre- 
sponding eigenstate is rather cumbersome. It will be sketched in appendix B. 



5.2.3 KPZ universality class 

From section 5.2.1, one immediately sees that the set of spectrum for the sector t (5.35) 
includes those for sectors consisting of single particles: 

Spec{si) C Spec(t), Si = {rm + ■ ■ ■ + mi} ^ ^m,+-+mi2mi+i+-+mN {1<1<N-1). 

(5.44) 

As discussed in section 3, the relaxation spectrum characterizing the universality class are 
the eigenvalues in the sector S;, whose real parts have the second largest value. As described 
below or in section 3, these eigenvalues form a complex-conjugate pair. We denote them by 

hereafter. The Bethe equation (5.33) describing Spec(s;) reduces to 



fe=l 



where we set the nesting order as (ai, 02) = (1, 2). Since the spectrum of the Hamiltonian is 
invariant under the change of the nesting order, and under the transformation p g^, it is 
enough to consider the case p > q and rii < L/2 



^This can be ecisily seen from the fact that the Bethe equation (5.33) is invariant under the transformation 
p <-> 5 and (ai, . . . , ojv) <-> (iV + 1 — oi, . . . , AT + 1 — oat). 
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For p ^ q, characterizes the KPZ universaUty class. The corresponding solutions to 
(5.45) are determined as follows [K]. Changing the variable A as 

pXj = , ^ = e-^ (5.46) 

we modify (5.45) and (5.34): 

l-e-4J =(-l)"-^n^',_e-2.^/ ^ = ^P-^^2^^[T^^- i-e-^ixj- 

(5.47) 

The meaning of 77 in the above will be revealed in section 5.3. Taking logarithm of both sides, 
one has 

± log ^/-^^^ ^ =1 ±y( logxk - log l-e-^^^fe/^i A (5_4g) 

where p = ni/L and /j G Z H ^-^^^ — . In fact, for sufficiently large ni and L, the following 

choice 



gives the solution corresponding to 



ni+l 



I o h ? for 1 < i < m — 1 , , , 



By carefully taking into account finite size corrections, the asymptotic form of E^ for 
L S> 1 is determined as 



E^ = ±2|(p -q){l- 2/>)KiL-i - 2C\p - p{l - p)L—2 + 0(L-2), (5.50) 

where C = 6.50918933794 ... [K]. Thus we conclude that the system for p / g belongs to the 
KPZ universality class whose dynamical exponent is 2; = 3/2. 



5.2.4 EW universality class 

For p = q, the set of eigenvalues of the Hamiltonian for an arbitrary sector t contains the 
relaxation spectrum corresponding to the "one-magnon" states. This can be seen by setting 
all the roots in (5.45) except for Ai to 1/p. Thus E^ are given by the second largest eigenvalues 
for this one-magnon states, and obviously do not depend on I. The Bethe ansatz equation 
determining the unknown Ai simply reduces to 

(pAi)^ = 1. (5.51) 

Solving this and substituting the solutions 

/ 27r/ci\ L 
pAi = exp(±-^j l<k<-. (5.52) 

into (5.34), we have E = —Apsm^{2-jTk/L). Obviously the case k = 1 gives the second largest 
eigenvalues: 

E^ = -4psin2 (^1) = -Apn'^L-'^ + 0(L-^), (5.53) 
which gives the EW exponent z = 2. 
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5.3 Pcirameterization with difference property 

Here we present the Bethe ansatz results in a more conventional parameterization [Sc, PS] 
with the spectral parameter having a difference property. 

First we treat the one-species case (N = 2) whose spectrum is given by (5.34) via the 
Bethe ansatz (5.45), where < ni < L and the nesting order is (01,02) = (1,2). Changing 
the variables as 

pXf'^ ^ exp{ipj + ??), ^ = e-^", (5.54) 



P 



we transform (5.34) and (5.45) to 

ni 

E = 2^^{cospj - A), 



e'^Pi = (-i)"i-ie-^^ n ^ , ^ , — ^=^, A = coshr/. (5.55) 



k=l 



This is nothing but the eigenvalue of the Hamiltonian for the XXZ chain threaded by a 
"magnetic flux" — ir/L: 



H = VpqJ2 je^^fe + + ^(cT^a^+i - 1)| . 



(5.56) 



The variable pj in (5.55) is called the quasi-momentum of the Bethe wave function. Intro- 
ducing the transformation (see [T] for example) 

1 Sh 

^^"^ ^=T^°^shi(u + 2)' ^^■=^K -1)' 

A(n) := - exp(ip(tx) + 77), ^ = A(iuf ^ - 1), (5.57) 

p J J 

we rewrite (5.55) in terms of the "rapidities" u^^^: 

ch(?7-u^ ') - ch?7 gi(«^ - 2i) 

where the two functions (f){u) and qi{u) are defined by 

. *w=n-'|(«-«f)- i^-'O) 

Applying the momentum-rapidity transformation (5.57) to the i2-matrix (5.2) (we write 
R{\{u)) = R{u)), we find the non-zero elements of R{u) can be written as 



RZiu) = 1, R%{u) = <; 



sh2fa+2) for a < /3, 



\T,tl for a > /?, 
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where a, /3 G {1, 2} in the present case. Up to the asymmetric factors e^^ and e^^"/^, these 
are the Boltzmann weights for the well-known six vertex model [Ba] associated with the 
quantum group Uq{sl{2)). For the gauge factors, see [PS, OY]. The i?-matrix (5.60) satisfies 
the Yang-Baxter equation 

i?23 (142)^13(^^1)^12(^^1 - U2) = Ruim - U2)Rl3{ui)R23{u2), (5.61) 

which possesses the difference property. The Hamiltonian (5.56) is expressed as the logarith- 
mic derivative of the transfer matrix T{u) (cf. (5.7)): 

T{u) = trwolRoLiiu -!)••• i?oi(i« - 1)], 



(5.62) 



2i^shr? d ~ 

H = — Inr('u) 

ri ou 

where N = 2 in the present case. Noting that 

d{\{m - 1)) = {pX{m - 1))^ = e"-^ (/>(«), 

1 sin -(u — V — 2i) 
f(X(m - 1), X(iv - 1)) = — — — — — = e"'' (5.63) 

the eigenvalue of the T{u) for the nesting order (01,02) = (1, 2) is given by 

Qi{u) qiiu) 

via the Bethe equation (5.58). 

The extension to the general (N — l)-species case is straightforward. We just let the local 
states a,P in (5.60) range over a,/3 G {1, • . . Finally, we write down the explicit form 

of the eigenvalues for an arbitrary nesting order: 

N—2 

+ (t)(u) V gfc(^-2i) gfc+l(^ + 2i) ^^L_^(n^^+^^_^^+2y:fe_,(«,_i-n,)e,fc+i-2nfc+i) 

+ . gAr-i(n-2i) ^,^_,(„^_^+2Efj-,^K-i-n,)g,._ (565) 
qN-i{u) 

Correspondingly the Bethe equation (5.33) is transformed to 

g-2,,(L-ni)ei2 ,/ (1)n ^ ^-r,L+r,(n2+2(ni-«2)) ^^^^j' ^ + ^^^^j' ^) 

' q^{uf-2i)q2{uf+2i) 
Q-'^v{J2j=i{nj-i-nj)eji+i-Y^'-jl\{nj-i-nj)dji) 

, ,0.- - ^^g^(4'^ + 2i)^^-l(^i•'^-2i) gz+i(uf) 

qi{uf-2\) qi-i{uf) g,+i(nf+2i) 

Q-'^v(.J2fji^{nj-l-nj)ejN-T,fS{'{nj-i-nj)ejN-l) 

, , qN-iiuf^~'^ + 2i) qN-2{uf-''^ - 2i) 

qN-i{u) - 2i) qN-2{u) ') 
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The spectrum of the Hamiltonian H is then determined by 



ch[r]Uj ) — chrj 
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A Proof of theorem 4.9 
A.l Mobius inversion 

The power set 5 (4.2) is equipped with the natural poset structure whose partial order is 
just C. In this appendix the partial order in M (2.15) induced via (4.3) will be denoted by 

^. Thus one has (4) ^ (1, 3) ^ (1, 2, 1) ^ (1, 1, 1, 1), etc. for L = 4. The description of < in 
Ai is pretty simple. In fact, those m' satisfying m' :< m = (mi, . . . ,m„) are obtained from 
m by successive contractions 

(. . .,mj,mj+i, . . .) H^- (. . . ,mj + m^+i, . . .). (A.l) 

Let ( = [Ci^' J s)) g/ be the \S\ x \S\ matrix defined by 

a5',5) = r (A.2) 

I otherwise. 

Since C is a triangular matrix whose diagonal elements arc all 1, it has the inverse /i = 
(/u(s', s)) ^, ^g^. fj, is called the Mobius function of 5, and is again a triangular (i.e., /x(s', s) = 
unless s' C s) integer matrix. 

Suppose /, y : (S — C are the functions on S. By the definition, the two relations 

/(s) = ^^5(5'), 5(5) = 5^Ms',5)/(V) iseS) (A.3) 

s'Cs s'Cs 

are equivalent, where the latter is the Mobius inversion formula. In a matrix notation, they 
are just f = gC and g = f/j,. In particular the sum involving //(s',s) can be restricted to 
s' C s. The Mobius function contains all the information on the poset structure. In our case 
of the power set S, it is a classical result (the inclusion-exclusion principle) that 

M(5^s) = (-l)t^'-«^ (A.4) 

where Jjs denotes the cardinality of s. 
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The Mobius inversion formula (A.3) and (A.4) on S can be translated into those on M. 
via the bijective correspondence (4.3). The result reads as follows: 



f(m) = ^ gim') (m £ M), (A.5) 
£?(mi,...,m„) = ^(-ir-V(ii, iimi,...,mn) eM), (A.6) 



where the sum in (A.6) extends over (ii, . . . , i;) G Ai such that (ii, . . . , i;) ^ (mi, . . . , m^). 
(We have written g{m) with m = (mi, . . . , m„) as g{mi, . . . , m^) rather than g{{mi, . . . , m„)), 
and similarly for /.) 

For m = (mi, . . . , m„) G corresponding to s G 5, we let m denote the clement in Ai 
that corresponds to the complement s = \ s G 5. Thus for L = 4, ~ acts as the involution 

0^ {1,2,3}, {1}^{2,3}, {2} ^{1,3}, {3} ^{1,2} 

on S, and similarly 

(4)^(1,1,1,1), (1,3)^(2,1,1), (2,2)^(1,2,1), (3,1)^(1,1,2) 

on A^. It is an easy exercise to check 



(mi, . . . , mn) = (i'"i-i2l'"2-221"*3-2 _ _ _ 21"^"-i-221'""-i) G M, (A.7) 
where "al~^6" should be understood as o + 6 — 1. 

A.2 Theorem 4.9 

We keep assuming the one to one correspondence (4.3) of the labels m £ A4 and $ E S and 
use the former. In view of (4.7) we have dim 14* = /(m) by the choice: 

/(m) = ( ^ ] := for m= (mi,...,m„) G X. (A.8) 

\^mi,...,m„y mi!---m„! 

Denote the g{m) determined from this and (A.6) by 

\mi, ...,mn/ 

Namely, we define 

(™,..':,™.>=E(-l)"-'(i„.^.,i,) for eA4, (A.10) 

where the sum runs over (ii, ■ ■ ■ E M such that (ii, . . . :< (mi, . . . , m^). From (4.7), 
(4.27), (A.5) and (A.8), we find dimX* = g{m). The function g{m) has the invariance 

mi,...,mn/ \mn,...,mi 

but it is not symmetric under general permutations of mi, . . . , m^ in contrast to /(m). 
Now theorem 4.9 is translated into 
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Theorem A.l. 

g{m) = g{rn) for any m G Ai. 

Example A.2. Take m = (L) hence m = (1-^). Then g{m) = = j = 1. On the 

other hand, g{fn) with L = 3 and 4 are calculated as 



\2,1 J \1,2 J \3 



^i,i,i,iy vAA) Vi'2,iy Vi'i'^y \3,\) \2,2j \i,3j \4 

Example A. 3. Take L = 5 and m = (1, 2, 1, 1) hence m = (2, 3). Then one has 
l,2,l,l) = (l,2,l,l) " (l,3,l) " (1,2,2) + (1,4 



3,1,1 J \3,2 J \4:,1 J \5 J ' 
2,3/ = (2,3) "(5 



(A.11) 
(A.12) 



The both of these sums yield 9. 
A.3 Proof 

We first generalize example A. 2 to 
Lemma A. 4. 

1)^(^1^) ^'^y L > 1. 

Proof. The left hand side is 1. From (A. 10), the right hand side is given by 

where the latter sum extends over ii, . . . ,ii G Z>i such that ii + ■ ■ ■ + ii = L. Thus we have 
the following evaluation of the generating functions: 

E%^^--M)'E '-i;;:::!r'" -ft-e-')'. 

L>1 

= E E = Ed - <'-^)' = o' - 1- 

L>1 ■ 1=1 1>1 L>1 ' 1>1 

The last relation tells that Y4=i Al^i = 1 ioY L >1. □ 
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1,4; \2,i,i 



In (A. 11), note that the first line of the right hand side equals ( , "^i ) ( o f i )' whereas 
the second line is nothing but ~ ( 3 ^ ]^ )' therefore one has 



1,2,1,1) (1,4) (2,1,1) (3,1,1 
The following lemma shows that such a decomposition holds generally. 
Lemma A. 5. 

L \ I L \ / L — nil \ / L 



mi,m2,. . . ,mn / \mij \m2, ■ ■ ■ ,mn / \mi + m2, ms, . . . , m„ / 

where ( ^ ) = ( ^ ,„. ) denotes the binomial coefficient. 

\rni/ \rni,L—mi/ 

Proof. In making (ii, . . . , ii) from (mi, . . . , m„) by the successive contractions (A.l), classify 
the summands in (A. 10) according to whether mi has been contracted to m2 or not. If 
it is not contracted, then ii = mi always holds and the corresponding summands yield 

L mi \ ^j^^ other summands correspond to the contracted case ii > mi + m2, 
mi/ \m2,...,mn/ 

whose contribution is — ( \ . □ 

\mi + m2,m3, . . . ,m„ 



Lemma A. 6. 



L 



where tt is an arbitrary array of positive integers summing up to L — t. 
Proof. We employ the induction on a. The case a= 1, i.e., 



follows from lemma A.5. Suppose that (A. 14) holds for a = a. Substitution of (A. 13) gives 

(l*,7r) "^^"^^'^"(s) _(a + l-s) (l*-«-i,7r ) ~ (a + 2-s,l*-«-2,7r)_ " 

(A.15) 

In the first term, the s-sum is taken as 

Thus (A.15) implies the a ^ a+1 case of (A. 14). □ 
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+ \ ^a-l 



Lemma A. 7. 

L 

where tt is an arbitrary array of positive integers summing up to L — a — 1. 
Proof. The case t = a + 1 in (A. 14) reads 



L 



L — a 
l,7r 



(l«+i,7r)"^^ (s) (a+l-s,7^)• 



Similarly, by setting t— s-a — l,a— ^a — 2 and tt — > (2, vr) in (A. 14), we get 

a-2 



a-2 



s,2,7r 



s=0 



L- s 
a — 1 — 8 



L-a + l 
2,7r 



a + 1 — S, TT 



where (A. 13) has been apphed in the second equahty. The sum of the two expressions gives 



^a+l ^ ) + \ ia-l 



L 

r-\2,7r 



a- 1 



L-a + l 



2,vr 



+ 



+ 



a-2 ^ 
s=0 ^ 



a — 1 — s 



L — a 
l,7r 

L-a+l 
2,7r 



The s-sum in the last term is {^]) due to (A. 16), finishing the proof. 
Proof of Theorem A.l. For (mi, . . . , m„) G A^, we are to show 



□ 



L 

mi, . . . ,m^ 



-|^mi-l2]^m2-2 2m„_i-22]^m„-l 



where the exphcit form of the dual is taken from (A. 7). We invoke the double induction on 
(L, n). The case L = 1 is trivially true. In addition, the case (mi, . . . , m„) = (L) has already 
been verified for all L in lemma A. 4. We assume that the assertion for {L',n') is true for 
L' < L with any n' and also for {L,n') with n' < n. Consider the decomposition (A. 13). By 



the induction assumption, the two quantities ^ y on the right hand side can be replaced 
with their dual. Then the relation to be proved becomes 
L 

"l^mi— l2]^"i2— 2 -|^m„-i— 22]^m„— 1 

L 
mi 



L — mi 

-|^m2 — l2]^"i3— 2 jm„-i— 22]^mri — 1 



2mi+m2 — l2J"^3— 2 jm„_i— 22jm„ — 1 

In terms of tt = (i'"2-221"^3-2 ;Lm„_i-22im„-i)^ ^j^jg jg expressed as 



L 

l"*i-\2,7r 
The proof is finished by lemma A.7. 



L 
mi 



L — mi 
l,7r 



□ 
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B Derivation of the stationary state 

Here we sketch a procedure to derive the stationary state in the framework of the algebraic 
Bethe ansatz. Throughout this appendix, we set the nesting order as a standard one: aj = 
j{l<j<N). 

As seen in section 5.2.2, the eigenvalue of the transfer matrix for the stationary state can 
be simply calculated by setting all the Bethe roots equal to 1/p. In contrast to the eigenvalue 
problem, the evaluation of the eigenstate is not trivial. This is caused by the S-operators 
such as (5.12) and (5.27) that approach zero as A — 1/p. Thus to obtain the state, we must 
normalize S-operators as 

i?S(A):=3^, a,G{j + 2,...,iV} (1 < < < j < iV - 2). (B.l) 

Note that B(°\\) := B{X). First we consider the eigenstate |{A(^-i)}) of r(^-2)(A|{A(^-2)}) 
As shown in section 5.1, this state is constructed by a multiple action of B^~^^ on |vac^^~^^) = 

|Ar-i)i®...® |Ar-i)„^_^: 

|{A(^-)})= ^ lim <-^)(Ar^^)---<-^^(Ai^:.^^)|vac(^-^)) 



"JV-i 



E "I n M (8) I"- 1+ E 

j=l ~ P^^] I V k=l j=l 



l<"/l< — <"/nff_l<nN- 

=: F(^-')"i-""^-2|ai)i ® • • • ® \an^_,)n^_„ aj e {N - 1,N} {1 < j < nN-2), 

(B.2) 

where element of |{A(^"^)}), from which the eigenstate KA^''^"^)}) of 

r(^-3)(A|{A(^-3)}) can be constructed as 

|{A(^-^)}) = hm f(^-3)"i--.-.s(^-3)(a(^-^)) . . . 5(^-3) (A(^:^))|vac(^-3)). 

a(^-2)^i/p 



(B.3) 



Note that denotes defined by (5.17). Since the coefficient f(^~^)"^ "''"Jv-2 

contains the term YI{1 —pXj^~'^^)~^, we cannot take the limit A^-^~^^ —>■ 1/p {1 < j < njv-2) 
independently of j. To take this limit correctly, we solve for the roots A^-^"^"* {2 < j < njv-2) 
in terms of A^^~^^ by using the Bethe equations (5.33): 

I = - TT - for iV > 3 (B 4) 

lljt=i JK'^j ^\ ) k=i J^^j ^\ ) 
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In the above, we have put A^- 



^) — I To extract the behavior of A^-^ around the point 
we expand them in terms of A^'^ — 1/p as 

njv-2 



fe=i 



p 



l/p)^ + 0{{X[ 



(JV-2) 



(B.5) 



where g[''^ = 6ik- Inserting them into (B.4) and comparing the coefficients of each order, one 
obtains the set of equations determining the coefficients gj'^^ . In the following, as an example, 
we write down the equations determining the first three coefficients. 



5(1) 



"•JV-2 



.(2) 1 



k=l %■ k=l { (^jk ^ 




l<fc<;<njv-2 I S'fc 



njv-2 



fe=i 



where Cj^ and 6 are defined as 



(1) 



(1) 



P9 



(1) 



r J3) 



J2) (2) 



(2) 



.(3) 



(2) , (1) (1) 

P9j +Pq9j 9k 



<^=,9f^-P9r+P,i9f^9?^+9P9^^ 



and, for N = 3, 



6^ = 1, b^'^=Lpgf, b('^=Lpgf+'-LiL-l)ipgff, 



and, for AT > 3, 



1, ^^^^ = E 



J,(3) 



fc=i 

1 - gAf 1 - gA;^-^) 



E 



fe=l U-P^fe 



(l-pAr-))2 ^^^^^ 



(B.6) 



(B.7) 



(B.8) 



(B.9) 



By solving these equations, the coefficients gj''^ are uniquely determined. For instance, gj^^ 
is simply given by roots of unity: 



exp 



27ri 
nN~2 



(j-l)} (l<j<niV-2). 



(B.IO) 
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Substituting (B.5) together with the explicit form of the coefficients g^p into (B.3) and then 

taking the hmit A^^~^^ — 1/p, one obtains the eigenstate |{A^^~^^}) whose elements give 

^(Ar-4)ai---a„jy_3 ^ Repeating this procedure, one calculates the stationary state. 

As a simple example, let us demonstrate this procedure for the maximal sector of L = 
N = Z{nj = L-j (1 < j < 2)). From (B.2) lA^^)) is given by 

1^^'^) = -P H) ® 1^)2 -P ^|3)i ® |2)2. (B.ll) 

Namely and F^"^ are, respectively, given by 

1 - pA^ ' l-pX\' 
To take the limit in (B.3), we expand Ag^' by solving (B.6). The resultant expression reads 

A^^^ = \ - {X? - I/P) + ^fcL^(AW - 1/pf + OiiX? - l/pf). (B.13) 

Inserting this and (B.12) into (B.2), and taking the limit A^^^^ — 1/p, we finally arrive at 
|{A(i)}) = -pip + g)C°i°2'^3|^^)^ ^ 1^2)2 ® 103)3, 

^123 ^ ^231 ^ ^312 = 2p + q, C^^^ = C^^^ = C^^^ =p+2q. (B.14) 

We leave it as a future task to extend the concrete calculation as above to the general 
case and derive an explicit formula for the stationary state. For an alternative approach by 
the matrix product ansatz, see [PEM]. 



C Eigenvalues and Bethe roots for N = L = 4 

We list the spectrum of the Hamiltonian (2.6) and the transfer matrix (5.5) in the basic 
sectors for iV = L = 4 and {p,q) = (2/3,1/3). The corresponding Bethe roots with the 
standard nesting order are also listed here. The spectrum of Hamiltonian is also obtained 
by specializing the result in figure 5. The second and the third examples from the last 
demonstrate that there are two Bethe roots that yield the same eigen-polynomial. These 
Bethe roots are not regular and the eigen-polynomial is not genuine in the sense of section 
5.1.2. The sectors are specified by (mi, . . . , 7714) according to (2.14). 
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Table 1: Table of eigenvalues of the transfer matrix and the Hamiltonian and the correspond- 



ing Bethe roots 1. 







{aW} 




{A(3)} 


(4,0,0,0) + i 














(1,3,0,0) M^ + i 





{1.5, 1.5, 1.5} 








lOX* 1 2A^ 1 ox 1 
-^27- + — 3- + - i 


-2 


( -2.51978 1 ft, ft, 
) 1.00989 - 0.,S6,'')2fi,'')i S V V 
[ 1.00989 + 0.,S6,'')2fi,'')i J 


(If ^ It) ^ ^ 5) 

+ (3Ti)^' + (3±0ATi 




-0.809148 ± 2.91994i ") 
0.824307 ± 0. 1825291 V 
1.16517 =F 0.618862i J 








(2,2,0,0) ^ + 1 





{1.5,1.5} 








/'32 1 4i N a4 2iA^ , 5A^ 

±iA=Fi 


(-If { 


-0.2411,58 ± 1. 941731 1 (?)(?) 
1.14116 ^ 0. 1417291 J V V 


4Al _ 2Ai ^ 3;^2 _ 3;^ + 1 


-3 


1 -0.75 - 1. 299041 \ 
\ -0.75 + 1.29904i J 








28A* 1 8A'^ 2A^ , 4A -i 
81 ^ 27 9 ' 3 


4 
3 


/ -3.62132 1 
\ 0.62132 / 








28A'' _|_ lOA^ 5A^ 1 5A i 
81 27 9 3 


5 
3 


f -2.42705 1 
\ 0.927051 / 








(3,1,0,0) ^ + 1 





{1.5} 








8A* , 8A^ 4A^ , r,> -1 


-2 


{-1.5} 








/'32_8i\i,4,('4 i4i\->3 

+ (iT|)A^-(|±OA±i 




{=Fl.5i} 


It) 


V) 


(1,1,2,0) ^+1 





{1.5, 1.5, 1.5} 


{1.5, 1.5} 





2A^ , 2A-' 2A^ .o\ 1 
— + -g 3-+^A-l 


-2 


C 1.00989 - 0. 5652651 I 
< 1.00989 + 0.565265i > 
[ -2.51978 J 


{1.5, 1.5} 


(71 


(§1,='= w) ^"^ ~ i4f ^ i) 
+ (l=FA)A2 + (|±i)A=Fi 




-0.809148 ± 2.91994i I 
0.824307 ± 0.182529i \ 
1.16517 =F 0.618862i J 


{1.5, 1.5} 


(71 

V 


('2 1 4i N 1,4 -r 2iA^ , 5A^ 

±iA=Fi 


(-1)^ 1 

I 


-0.241158 ± 1.94173i ] 
1.14116 =F 0. 1417291 > 
1.5 J 


f 0.69207 ± 0. 6483161 \ 
\ 1.19884 ^ 0. 1612861 J 





- ^ + 3A2 - 3A + 1 

81 3 


-3 1 
I 


1.5, -0.75 - 1.29904i 1 
-0.75 + 1.29904i J 


/ 0.789474 - 0.957186i \ 
\ 0.789474 + 0.9571861 / 





14A'' , 8A^ 2A^ 1 4A i 

81 ^ 27 9 3 


4 
3 


f 1.5, -3.62132 1 
\ 0.62132 J 


/ 0.686441 - 0. 5033641 1 

\ 0.686441 + 0. 5033641 J 





14A'* I lOA^ 5A^ 1 5A -, 
81 27 9 ' 3 


5 
3 


f 1.5, -2.42705 1 
\ 0.927051 / 


f 0.765957 - 0.1888111 1 
\ 0.765957 + 0.1888111 J 





_ 22d + 4Ai _ 2A - 1 


1 

' \ 


f -0.61352 - 2.()(i5:i(ii 1 
[ ().727U4 J 


f -9.97723 1 
I 0.977226 / 





(^±M)A^ + (iTi)A3 
+ (i±f )A2-(i±3i)A±i 


-3±i { 


-1.97293 =F 0.68627i ") 
-0.405293 ± 1.74466i \ 
0.762834 T 0. 4814631 J 


/ -3.45362 ± 7. 737281 1 
\ 1.07431 T 0. 2890051 / 





(1,2,1,0) ^ + 1 





{1.5,1.5,1.5} 


{1.5} 





~81~ + ~9 3~ + ^A - 1 


-2 


f 1.00989 - 0. 5652651 I 
I 1.00989 +0. 5652651 ) 
[ -2.51978 J 


{1.5} 





(^,='= It) ^ =^ i) 

+ (lTi)A2 + (|±i)ATi 


{ 


-0.809148 ± 2.91994i "I 
0.824307 ± 0. 1825291 [. 
1.16517 ^ 0. 6188621 J 


{1.5} 





lOA'' 1 8A-' 4A^ 1 r,\ -, 


-2 


(-1.5, 1.5, 1.5} 


{-0.5} 





(|Tf)A4 + (^±f)A3 
+ (§T|)A^-(l±OA±i 




{=Fl.5i, 1.5, 1.5} 


{0.253846 =F 0.9692311} 






42 



Table 2: Table of eigenvalues of the transfer matrix and the Hamiltonian and the correspond- 



ing Bethe roots II. 


sector A{X) 


E 


{AW} 


{Ai2)} 




(1,2,1,0) + ^ + 


-1 


f -7.5494 1 
^ 0.383413 ) 
{ 1.16599 J 


{1.5} 





- ± 3i) A2 ± 3iA =F i 


(-3)2 


f -1.9214 ± 0. 2385121 1 
i -0.315515 ^ 1.62413i !■ 
[ 1.00162 ± 0. 3267981 J 


{-2.19231 =F 3.461541} 





4A 1 'A 1 14A 'All 
27 27 9 3 "f" 


7 

~3 


f -0.511746 - 1.74231i ") 
^ -0.611746 + 1.74231i ) 
( 1.02349 J 


{1.5} 





2A* 4A^ _|_ 20a2 8A i , 
9 27 9 3 


8 
3 


f -0.717003 - 1.5691i I 
< -0.717003 + 1.5691i ) 
1 1 1 34m 1 


{-0.214286} 





(2,1,1,0) a|Al + i 





{1.5, 1.5} 


{1.5} 





/ 8 1 4i \ x4 -,- 2iA3 , 5A^ 
1.27 =■= 8TJ T -g- + -g- 

itiA ^ i 


(-1)^ 


f -0.241158 ± 1. 941731 \ 
1 1.14116 =F 0. 1417291 J 


{1.5} 





+ ■■-lA- - :U - 1 


-3 


/ 0.75 1.2<)<)().lj \ 
\ -0,7.'; -r i. 299041 J 


{1.,-,} 





20A* I 8A^ 2A^ , 4A i 
81 27 9 3 


4 
3 


f -3.62132 1 
\ 0.62132 / 


{1.5} 





20A* , lOA^ 5A^ , BA i 
81 27 9 3 


5 
3 


r -2.42705 1 
\ 0.927051 / 


{1.5} 





4A'' , SA^ 4a2 , 1 
— + -g — +ZA- 1 


-2 


{-1.5,1.5} 


{0.3} 





(iTt)A4 + (A±|)A3 
+ (§=F|)^'-(i±i)A±i 


1 J- i 


{=Fl.5i, 1.5} 


{0.617647 =F 0. 5294120 


It) 


^+8Ai_^2Ai_2A+l 

81 9 3 


-2 


{-1.5i, 1.5i} 


{-4.5} 





_(|± 2ft)A2 + (i±3i)Aq=i 


-3±i 


{-1.5,=Fl-5i} 


{-2.55882 =F 4.764711} 





(1,1,1,1) ifl + i 





{1.5,1.5,1.5} 


{1.5,1.5} 


{1.5} 


10Al+2Ai_2Ai+2A-l 


-2 


C 1.00989 - 0.5652651 } 
I 1.00989 + 0.5652651 1 
-2.51978 J 


{1.5,1.5} 


{1.5} 


^ It) ~ ^ i) -''^ 

+ (§=Fi)A2 + (i±i}ATi 




C -0.809148 ± 2.919941 } 
■! 0.824307 ± 0.1825291 [• 
[ 1.16517 =F 0.6188621 J 


{1.5,1.5} 


{1.5} 


/lO 1 4i \ >4 2iA^ , 5A^ 
l8T='=8T-''^ T~9~ + ~9~ 
±iAq:i 


(-1)2 


f -0.241158 ± 1.941731 ") 
^ 1.14116 =F 0.1417291 > 


r 0.69207 ± 0.6483161 1 
\ 1.19884 =p 0.1612861 J 


{1.5} 


^ - ^ + 3A2 - 3A + 1 


-3 


C -0.75 - 1.299041 "1 
■j -0.75 + 1.299041 | 


\ 0.789474 + 0.9571861 J 


{1.5} 


2A* , SA'"* 2A^ 1 4A -, 
27 27 9 ' 3 


4 
3 


/ 1.5, -3.62132 1 
\ 0.62132 / 


/ 0.686441 - 0.5033641 \ 
\ 0.686441 + 0.5033641 J 




2A* , lOA^ 5A^ , 5A i 
27 27 9 3 


5 
3 


/ 1.5, -2.42705 1 
\ 0.927051 J 


r 0.765957 - 0.1888111 1 
\ 0.765957 + 0.1888111 / 


{1.5} 


^ + ^ + ^-2A + l 


-2 


C -0.61352 - 2.065361 ") 
I -0.61352 + 2.065361 ) 
[ 0.72704 J 


/ -9.97723 \ 
\ 0.977226 / 


{1.5} 


+ (i±fi)A2-(i±3i)A±i 


-3±l 


f -1.97293 =F 0.686271 } 
{ -0.405293 ± 1.744661 \ 
[ 0.762834 =F 0.4814631 J 


f -3.45362 ± 7.737281 1 
\ 1.07431 T 0.2890051 J 


{1.5} 


# + ^ + A-l 


-1 


/ 1.16599, -7.5494 1 
\ 0.383413 J 


{0.5, 1.5} 


{0.954545} 


(-i?Tt)A^+(f 

-(§ ±3i) A2±3iA=Fi 


(-3)2 


f -1.9214 ± 0.2385121 ") 
i -0.31BB1B =F 1.624131 V 
( 1.00162 ± 0.3267981 J 


/ -2.19231 =F 3.461541 \ 
1 15/ 


{-0.101751 =F 0.590811} 
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Table 3: Table of eigenvalues of the transfer matrix and the Hamiltonian and the correspond- 



ing Bethe roots III. 



sector 


A(A) 




{Ail)} 






{Ai^^} 


(1,1,1,1) 

1j2 L 

27 ^ 


14A^ 7A 1 1 
_ 3" + J- 


7 

3 


f -0.611746 - 1. 742311 1 
I -0.611746 + 1.74231i ^ 
( 1.02349 J 


{-1.5,1.5} 




{0.3} 


4A _|_ 

27 27 


ZUA oA 1 1 

9 3 


8 


f -0.717003 - 1.5691i "I 
I -0.717003 + 1.5691i } 
( 1.13401 J 


J — [J.ZL^ZOO 1 

1 15/ 




{0.672414} 


81 ^ 9 ' 


f ^ - 2A + 1 


-2 


{-1.5i, 1.5i, 1.5} 


{-2.06427, 0.778553} 




{-1.77273} 




+ (44 ± I) x-- 




{-1.5.=1.5i.l.5| 


/ -2..3.")1.").") = :-i.()2941i 1 


{ l.il:'/ 


roil : l.')1878i} 


81 9 


_ 4Ai + 2A - 1 


-2 


{-1.5,1.5,1.5} 


{-0.5, 1.5} 




{0.576923} 








{1.5, 1.6} 




{-1.5} 


(^„'F |r) -"^^ 

+ (|t|)a^ 


+ (^±f)A^ 
-(|±i)A±i 




{=Fl.5i,1.5, 1.5} 


/ 0.253846 =F 0.969231i \ 
I 1-5 / 


{0.784404 =F 0.3853211} 










{1.5,1.5} 




{=Fl.6i} 




- 6A2 + 4A - 1 


-4 


{-1.5, -1.5i, 1.5i} 


f -0.5 - 3.74166i 1 
\ -0.6 + 3. 741661 / 




{25.5} 
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